Search for Neutral Q-balls in Super-Kamiokande II 
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A search for Q-balls has been carried out in Super-Kamiokande II with 541.7 days of live time. A 
neutral Q-ball passing through the detector can interact with nuclei to produce pions, generating 
a signal of successive contained pion events along a track. No candidate for successive contained 
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event groups has been found in Super-Kamiokande II, so we obtain upper limits on the flux of such 
Q-balls. 

PACS numbers: 95.35.+d, 14.80.-j 



In the framework of some supersymmctric models, sta- 
ble non-topological solitons called Q-balls [l[ can be pro- 
duced in the early universe and contribute to the dark 
matter 0, H, 13) B @| • If Q-balls are formed in the early 
universe and have survived until the present time, they 
would be part of the dark matter located in the halo of 
Galaxy. Q-balls are solitons with theoretically appealing 
roles in some dark matter scenarios and in explanations 
of the baryon asymmetry Q • 

In this model [3j, the Q-ball mass is given by 



and its radius is 
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where Ms is the energy scale of SUSY breaking and Q 
is the Q-balPs baryon number. Q-balls are classified into 
two groups according to the properties of their interac- 
tions with matter: supersymmetric electrically charged 
solitons (SECS) and supersymmetric electrically neutral 
solitons (SENS). Scintillator detectors may be more suit- 
able for the detection of charged Q-balls with large energy 
loss in matter. Thus, we will confine our attention here 
to neutral Q-balls. 

When a neutral Q-ball collides with a nucleon it ab- 
sorbs its baryonic charge and induces the dissociation 
of the nucleon into free quarks. In this process, about 
1 GeV energy is released by the emission of typically two 
or three pions [8, 9]. This is called the KKST process 

The cross section of the interaction between neutral 
Q-balls and the matter is roughly estimated to be the 
geometrical size of the Q-ball, 



er = ttRq 
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We consider cross sections between 0.02 and 200 mb in 
this work. Neutral Q-balls lose energy due to collisions 
with nucleons, and the rate of energy loss is Q, 
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where £ = 1 GeV is the released energy in one collision 
and A is the mean interaction length. For 100 mb cross 
section, the energy loss of Q-ball is estimated to be about 
60 MeV/cm. The mass of a Q-ball with this cross section 
is about 10 19 MeV/c 2 . Therefore the effect of the energy 
loss in the Earth is expected to be negligible compared 
to the kinetic energy of a Q-ball. 

Super-Kamiokande is a cylindrical 50kton water 
Cherenkov detector located at a depth of 2,700 m wa- 



ter equivalent. The water tank is optically sepa- 
rated into two regions; the inner detector (ID) instru- 
mented with inward facing 20 inch diameter photomu- 
tiplier tubes (PMT) and the outer detector (OD) in- 
strumented with outward facing 8 inch PMTs. Super- 
Kamiokande II (SK-II) is the second phase of Super- 
Kamiokande experiment and started taking data in De- 
cember 2002 with 5J82 PMTs in the ID and 1,885 PMTs 
in the OD. See [l5j for more details on the detector. 
The inner PMTs are instrumented with acrylic and FRP 
(fiber reinforced plastic) covers to avoid a chain reaction 
implosion. The event gate width is 1.3/isec. The OD is 
used to veto entering cosmic ray muons and to tag exit- 
ing charged particles. In this paper we report the results 
from a search for neutral Q-balls using 541.71ive days of 
data for the SK-II period. Previous searches for WIMP 
dark matter with SK-I data have been reported in [Icj |. 
A Q-ball interacting in the detector will produce several 
pions which emit Cherenkov light. Thus, we can observe 
these pion events along the Q-ball track in the detector. 

We developed a Monte Carlo simulation to estimate 
the detection efficiency of Q-ball events. The direction 
of a Q-ball is randomly generated, and the distance be- 
tween the center of the SK detector and a Q-ball track is 
assumed to be uniformly distributed within 50 m. The Q- 
ball velocity is generated assuming a Maxwellian distri- 
bution with a mean value of 270 km/sec — a typical galac- 
tic velocity. Therefore, the velocity of Q-balls observed 
on the earth is determined by superposing the velocity 
of solar system 220 km/sec on the Q-ball velocity in the 
halo [HI, [IH. The pions produced in the interactions 
with a Q-ball and a nucleon can be determined using 
KNO scaling [l3| . The total energy of the released pions 
is equal to the nucleon mass. The charges of the pions 
are randomly generated, so that the probability of ir + , 
7r~ , 7T° generation is equal within the constraint of the 
charge conservation. The kinematics of generated pions 
are determined by assuming the Q-ball and nucleus scat- 
ter into the Q-ball and the several pions. The number of 
generated Q-ball tracks is 20,000 for each cross section. 

In this work we look for the signal of a neutral Q-ball 
interacting successively with two or more nucleons. The 
main strategy of the data reduction is to reject muon 
events in order to observe a Q-ball signature consisting 
of at least two successive events in which 1 GeV of energy 
is released in the form of pions in the ID. A Q-ball with 
/3=10~ 3 covers 30m in 100 /isec; assuming cr=10mb it 
interacts on average 18 times. The details of the event 
reconstruction in Super-Kamiokande can be found in [lij . 

The first selection criteria are: 

1-1) The number of total photoelectrons (p.e.'s) in the 
ID should be larger than 300 and the maximum num- 
ber of p.e.'s in the ID within sliding 300nsec time win- 
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FIG. 1: The charge distributions for the data and Q-ball 
Monte Carlo events. The solid black line indicates the data 
and the dotted red line indicates the Monte Carlo events with 
cr=10 mb. 



dow should be less than 10,000, which corresponds to 
about 60MeV and 2GeV energy for electrons, respec- 
tively. These cuts reject electron events from decays 
of cosmic ray muons, high energy through-going muon 
events and electrical noise events. 

Figure [1] shows the charge distributions for the data 
and Q-ball Monte Carlo events with cr=10mb. The total 
number of p.e.'s in the ID is shown in the upper figure. 
The rising tail below a few hundred p.e.'s in the data 
has a lot of low energy events from radioisotopes and 
electrical noise events. On the other hand, for the Q- 
ball Monte Carlo events, the rising tail below 100 p.e.'s 
comes from the events which do not enter the ID. The 
peak around 200 p.e.'s consists of electrons from decay 
of charged pions. The events of neutral pions form the 
peaks around 4000 p.e.'s. In the bottom figure we show 
the number of ID p.e.'s within a 300nsec time window 
for events with more than 300p.e.'s in the ID.. 

1- 2) At least one delayed event should exist within a 
100/xsec window after the event selected by 1-1). This 
cut rejects single events. 

All events selected by 1-1) and 1-2) are considered to 
be an event group. The event selected by criterion 1-1) 
is regarded as the first event in an event group. 

After the first selection, there are still background 
event groups which are caused by cosmic ray muons, elec- 
trons from a decay of cosmic ray muons and electrical 
noise caused by large muon pulses. Background event 
groups still survive the above selection criteria with an 
approximate rate of 0.06 /sec. 

In the second reduction, we require at least two events 
which are entirely contained in the ID to reject cosmic 
ray muon events. At least two events in an event group 
must satisfy both of 2-1) and 2-2) criteria: 

2- 1) The number of hit PMTs in the largest OD hit 
cluster should be less than 16, in order to remove pri- 
marily cosmic ray muon events as shown in Figured) 




number of hit PMTs in the largeset OD hit cluster 



FIG. 2: The distribution of the number of hits in the largest 
OD cluster for the data and Q-ball Monte Carlo events. The 
lines are same as in Figure [1] 
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FIG. 3: The distribution of the goodness of the vertex fit for 
the data and the Q-ball Monte Carlo events. The lines are 
same as in Figure [T] 



2-2) The number of p.e.'s in the ID within 300nsec 
time window should be larger than 500 and the goodness 
of vertex fit should be larger than 0.6 (see Figure [3]) to 
better reject events associated with electrical noise. 

2-3) In a event group two or more events should satisfy 
both 2-1) and 2-2) criteria. 

Two event groups have been found, both containing 
two successive events. However, they are accidental co- 
incidence events consisting of a low energy event and a 
cosmic ray muon within the 1.3/xsec trigger gate. The 
cosmic ray muon is separately recorded over two events 
since the cosmic ray muon hits at the end of the trigger 
gate. These event groups are not rejected in the former 
reduction stages because of the absence of hit OD PMTs 
and the large number of total p.e.'s in the ID due to the 
muons. 

We design the following selection criteria to reject the 
accidental coincidence events without losing any Q-ball 
Monte Carlo events. 

3) If an event group contains two events, for the first 
event, the ratio of the number of hit ID PMTs in a 
lOOnsec time window at the end of the trigger gate to 
the total number of hit ID PMTs is required to be less 
than 50%. The distributions of this ratio for the data 
and the Monte Carlo events are shown in Figure Q] The 
events for which the ratio is more than 50 % are defined 
to be the accidental events. 
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FIG. 4: The distribution of the ratio of the number of hit ID 
PMTs in a lOOnsec time window at the end of the trigger 
gate to the total hit ID PMTs for the data and the Monte 
Carlo events. The lines are same as in Figure [T] 
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TABLE II: Summary of flux upper limits on neutral Q-balls 
for several Q-ball-nucleon interaction cross-sections. 



Reduction step Data Monte Carlo (10 mb) 



1-1) 
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2-1) 


2610998 


2865(14.3%) 


2-2) 
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TABLE I: Number of event groups after each reduction step 
for data and the Q-ball Monte Carlo with <r=10mb. 



No group of events has survived after applying all of 
the selection criteria. The number of background event 
groups due to the accidental coincidence of 2 or more 
atmospheric neutrino events within 100 usee during the 
541.7 day exposure is estimated to be 8.6xl0~ 5 . Table |l] 
shows the number of Q-ball candidates after each reduc- 
tion step for data and Monte Carlo with <r=10mb. 

Since there is no candidate event group, the upper lim- 
its on the Q-ball flux for various mean interaction lengths 
can be calculated from the following formula, 

Flux < (90%C.L.), (5) 

where A^ = 2.3 for no observed event group, T is the 
detector live time (541.7days), and S e ff£l, the effective 
aperture, which is defined as: 

S e fftt = enr 2 x 4vr [m 2 sr], (6) 

where e is the detection efficiency and r is the impact 
parameter, which is defined to be 50 m in this analysis, 
e can be estimated by Monte Carlo events separately for 
four different values of Q-ball cross section and is 4.5%, 
13.6%, 14.5% and 14.4% respectively for lmb, 10 mb, 
100 mb and 200 mb cross section. 

In another theory, Q-balls could convert matter into 
antimatter on their surfaces ||. In this baryon, 
reflected off a Q-ball as an anti-baryon, would quickly 
annihilate with a baryon and release about 2 GeV of mass 



energy. The present selection criteria is as efficient for the 
Q-balls from this theory. We also find that varying the 
multiplicity of the generated pions from two to four does 
not change the detection efficiency. The velocity of the 
neutral Q-ball has uncertainty of as much as ±10% [U, 
[l2l ]. The efficiency is not affected by this uncertainty. 
Table [TT1 summarizes the obtained flux limits. 

Figure \E\ shows the 90% C.L. upper limits on the neu- 
tral Q-ball flux as a function of the Q-ball cross sec- 
tion. The upper right hatched region is excluded be- 
cause the Q-ball mass density exceeds the Galactic mass 
density [l6|. The signal of the neutral Q-balls may be 
close to that of nucleon decays catalyzed by monopoles, 
so it is possible to obtain a rough estimate of the flux 
upper limit using the results of monopole search experi- 
ments .16]. The hatched region is previously excluded by 
the experiments which searched for monopole-catalyzed 
nucleon decays. The triangle points are the results from 
Kamiokande [l7| and the square point is obtained by 
MACRO The present results, shown in circles, are 
the most stringent limits on neutral Q-ball flux for Q- 
ball cross section below 200 mb. We require at least 2 
successive events, so the detection efficiency would be 
proportional to the square of the Q-ball cross section in 
the small cross section region. Assuming that the flux 
limits should be in inverse proportion to the square of 
cross section, the dashed line is obtained in the region 
below 1 mb cross section. 

In summary, we report on a search for groups of suc- 
cessive contained events at Super-Kamiokande II. We 
have found no evidence for the successive contained event 
groups in 541.7days of SK-II data. New upper limits at 
90% C.L. for neutral Q-ball flux are obtained. These 
limits are the most stringent bounds on Q-ball flux for 
Q-ball cross section below 200 mb. 

We gratefully acknowledge the cooperation of the 
Kamioka Mining and Smelting Company. The Super- 
Kamikande experiment has been built and operated from 
funding by the Japanese Ministry of Education, Culture, 
Sports, Science and Technology, the United States De- 
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FIG. 5: 90 % C.L. upper flux limits on the neutral Q-ball flux 
as a function of the Q-ball-nucleon interaction cross section. 
The circles show the results of this paper and the dashed line 
in the region below 1 mb cross section is obtained assuming 
that the flux limits should be in inverse proportion to the 
square of cross section. 
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